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Background: As one way of assessing the impact of Avahan, the India 
AIDS Initiative of the Bill & Melinda Gates Foundation, we examined the 
association between HIV prevention program indicators and changes in HTV 
prevalence among female sex workers (FSWs) between 2005 and 2009. 
Methods: We conducted a secondary data analysis from 2 large cross- 
sectional surveys (2005-2006 and 2008-2009) across 24 districts in 
south India (n = 1 1,000 per round). A random-effect multilevel logistic 
regression analysis was performed using HTV* as the outcome, with in- 
dividual independent variables (from both surveys) at level 1 and district-level 
FSW-specific program indicators and contextual variables at level 2. Program 
indicators included their 2006 value, the difference in their values between 
2008 and 2006, and the interaction between this difference and study round. 
Results: HTV prevalence among FSWs decreased from 17.0% to 14.2% (P 
< 0.001). This decline varied significantly (P = 0.006) across levels of dif- 
ference in program coverage (% of FSWs contacted by the program in a given 
year). Odds ratios comparing HIV prevalence between rounds changed with 
the level of increase in coverage and were statistically significant with 
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coverage increase > quartile (Q) 1: odds ratio, 0.85 at Ql; 0.78 at Q2; 0.66 at 
Q3; and 0.51 atQ4. 

Conclusions: These findings suggest that increased program coverage 
was associated with declining HIV prevalence among FSWs covered by 
the Avahan program. The triangulation of our results with those from 
other approaches used in evaluating Avahan suggests a major impact of 
this intervention on the HIV epidemic in southern India. 

A vahan, the India AIDS Initiative of the Bill & Melinda Gates 
f\ Foundation, is a comprehensive HIV prevention program 
targeting the most at-risk populations in the states of India most 
affected by the HIV epidemic (Andhra Pradesh, Karnataka, 
Maharashtra, and Tamil Nadu in the south of the country, and 
Manipur and Nagaland in the northeast). It was initially imple- 
mented through a network of nongovernmental organizations 
(NGOs) and has progressively been transferred to government 
agencies and community-based organizations since 2008. 

Because cluster randomized designs have major limita- 
tions for evaluating complex large-scale interventions, 1 ' 2 the initial 
impact evaluation plan for Avahan included serial cross-sectional 
studies in the high-risk populations, based on which sophisticated 
mathematical models were developed to estimate program im- 
pact in those populations, and in the general population. 1,3 The 
use of this modeling approach has been complicated by the lack 
of adequate baseline data, 4 but recently, published results using 
this approach suggest a major impact of Avahan on the HIV 
epidemic in Karnataka 5 ' 6 and in the 4 South India states men- 
tioned above. 7 A purely statistical approach based on HTV surveil- 
lance data among pregnant women was also used to estimate 
the overall Avahan impact in the 6 Avahan states. 8 Although 
suggesting a positive impact of the intervention in the general 
population, this study has the limitation of not taking into ac- 
count the nature of the Avahan intervention, including its var- 
ious components, and does not allow for the analysis of the 
actual mechanism leading to the impact. 9 

We examined changes over time occurring in HIV preva- 
lence among female sex workers (FSWs), the largest risk popu- 
lation in the program, between 2 large surveys carried out in 
2005 to 2006 and 2008 to 2009, in 24 districts of the 4 southern 
states of India covered by Avahan, where the HIV epidemic is 
predominantly heterosexual and mainly driven by sex work. 10 ' 11 
As indicators related to program coverage and implementation 
have shown very sharp increases between baseline and 2008, 12 
we considered the impact of these increases on the changes in 
HIV prevalence between surveys. Given the heterogeneity of 
the HIV epidemic in India, 13 it is important to also consider the 
influence of broader contextual factors that may have an impact 
on the HIV epidemic in high-risk groups. We therefore carried out a 
multilevel statistical analysis of the individual, contextual and pro- 
gram factors impacting on HIV prevalence among FSWs, and on 
the changes in the prevalence levels between the 2 survey rounds. 
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METHODS 

Data Sources 

The design and study procedures of the cross-sectional 
studies, known as integrated behavioral and biological assess- 
ments (IBBAs), have been described in detail elsewhere 14 ' 15 and 
are summarized in the Supplementary Online Material Methods 
section, http://links.lww.com/OLQ/A86, as are the main com- 
ponents of the Avahan intervention. 

The common set of program indicators for all districts 
covered by Avahan has been described extensively elsewhere. 12 
For the purpose of the present analysis, we used data collected 
by the NGOs implementing the intervention in each district 
along with complementary information to estimate: (1) program 
coverage (proportion of FSWs contacted by the program at least 
once in a given year), (2) the proportion of FSWs contacted 
monthly, and (3) the proportion of condom requirement met in 
a given year (number of condoms distributed in a given year 
divided by the estimated number of commercial sexual acts by 
all FSWs in the same year). Details on the methods used to 
derive these estimates are provided in the Supplementary Online 
Material Methods section, http://links.lww.com/OLQ/A86. The 3 
indicators were computed for years 2006 to 2008 because the in- 
formation was complete for all districts only for that period (the 
program started being handed over to the government in 2008). In 
the statistical models, we used the 2006 values of the indicators 
and the difference between the 2008 and 2006 values. 

We considered a total of 53 district-level contextual 
variables for inclusion in the multilevel statistical model (see 
complete list in the Supplementary Online Material Table, 
http://links.lww.com/OLQ/A87). 

The study and its consent procedures were approved by 
the ethics committees of all institutions that were involved in the 
data collection for this study (see Supplementary Online Mate- 
rial Methods section for more details on ethical considerations 
http://links.lww.com/OLQ/A86). 

Statistical Analyses 

We assessed differences in personal characteristics be- 
tween participants who were tested for HIV and those who 
were not, and between participants in round 1 and those in round 
2, using the x 2 test. We compared HIV prevalence between the 
study rounds for each of the 24 IBBA districts. We then built 
multilevel (FSW personal characteristics at level 1 and district- 
level contextual and program variables at level 2) logistic regres- 
sion models with HIV prevalence as the outcome, using the 
round 1 and round 2 merged data set. Study round was consid- 
ered as the main independent variable in our analyses because 
the difference in HIV prevalence between the rounds was our 
central interest. We also preselected a certain number of the 
individual variables from the questionnaire for inclusion as in- 
dependent variables in the model because they could be con- 
founders of the differences observed between rounds. This was 
done based on variables previously identified to be associated 
with HIV in the first IBBA round 14 and of a priori knowledge 
from the literature (see the Supplementary Online Table for the 
complete list of individual variables considered). We entered all 
of these variables into a multilevel logistic regression model 
with a random intercept for the districts and a random coeffi- 
cient for round because the difference in HIV prevalence varied 
widely between districts. We then removed sequentially from 
this model all the individual variables with P > 0.05, as long as 
they were not confounding the other associations, to arrive at a 
model with individual variables only. 

We preselected the contextual variables for addition to 
the models above based on linear regression models with the 



aggregate HIV prevalence rate in each district as the dependent 
variable and the contextual variables as the independent vari- 
ables. Those with a P value less than 0.05 in these models were 
then added to the previous multilevel models. Baseline HIV 
prevalence among FSWs in each district was also included in 
the model as an additional district-specific variable. The removal 
of nonsignificant variables at this stage led to the basic model with 
individual and district-level contextual variables. 

To the basic model, we then added the 3 pairs of district- 
level program variables (2006 value and difference between 
2008 and 2006), along with the interaction terms between the 
2008 to 2006 difference and the round. We kept in the model 
only the program variables corresponding to a significant in- 
teraction (P < 0.05) with the round. Finally, because we were 
interested in factors that could impact the differences in the out- 
comes observed between the IBBA rounds, we examined one- 
by-one the interaction terms between the round and each of 
the other individual and district-level contextual factors included 
in the first multilevel model described above. We kept the sta- 
tistical interactions that reached P < 0.05 in these models. This 
led to the final model, with assessment of the intervention impact, 
through the interaction terms between the study round and the 
2008 to 2006 differences in program variables, and inclusion of 
other factors having a multiplicative statistical interaction with the 
round. All analyses described above were carried out using 
STATA IC 1 1 . 1 for Windows (StataCorp LP, College Station, TX). 

RESULTS 

Among the 14,601 women approached in round 1, 1 1,690 
(80.1%) answered the questionnaire. Among the latter, 10,806 
were tested for HIV (92.4%), 10,716 (91.7%) for syphilis, and 
10,674 (91.3%) for gonorrhea and chlamydia (NG/CT). The 
overall participation rate was slightly lower in 2008: 15,682 women 
were approached, of whom 11,486 (73.2%) answered the ques- 
tionnaire; however, among the participants, 10,901 (95.0%) were 
tested for HIV, 10,806 (94.1%) for syphilis and 10,880 (94.8%) 
for NG/CT. No data were available to compare participants with 
nonparticipants. However, Table 1 shows the comparison be- 
tween participants who were tested for HIV and those who were 
interviewed only. There were significant differences for most 
variables examined. This is, however, largely due to the very 
large sample sizes of these surveys. In addition, the proportions 
among all participants were almost identical to those among 
those tested for HIV (see the "Total" columns in Table 1). 

Table 2 shows the comparison between various character- 
istics of participating FSWs between rounds 1 and 2. The differ- 
ences were significant for most variables examined. However, 
this was mostly for behavior or infections that were targeted by 
the Avahan program, such as condom use, violence, NG/CT, and 
syphilis prevalence. For the other significant variables, the actual 
differences in the distribution of their values between rounds 
were mostly small, with the statistical significance often due to 
the very large sample sizes of these surveys. 

HIV prevalence declined significantly from 17.0% to 
14.2% (P < 0.001) from 2005-2006 to 2008-2009 (Table 3). 
This decline was significant in 2 of the 4 states and 10 of the 
24 districts. There was an increase in HIV prevalence in 4 districts, 
but it was statistically significant in only 1 . Table 3 also shows 
the variation in program coverage and condom requirement met 
between the study rounds. There was a significant correlation 
between changes in HIV prevalence between 2006 and 2008 and 
changes in coverage (Spearman r = — 0.410, P = 0.048), whereas 
such an association was borderline significant when consider- 
ing condom requirement met (Spearman r = —0.398, P = 0.054). 
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TABLE 1 . Comparison of the Characteristics of Female Sex Workers Providing and Not Providing a Biological Sample Among Those 
Interviewed in Survey Rounds 1 and 2 
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Variable 


% 


% 


% 




% 


% 


% 




Current age, y 


















15-24 


24.7 


20.9 


21.2 


O.001 


22.6 


17.8 


18.0 


<0.001 


25-34 


52.0 


45.9 


46.3 




47.5 


44.3 


44.4 




35+ 


23.3 


33.3 


32.5 




29.9 


38.0 


37.6 




Literacy 


















Illiterate 


70.8 


65.1 


65.6 


0.001 


52.1 


60.7 


60.2 


<0.001 


Literate 


29.2 


34.9 


34.4 




47.9 


39.3 


39.8 




Source of income other than sex work 
















No 


76.2 


57.1 


58.5 


<0.001 


62.0 


54. 1 


54.5 


<0.001 


Yes 


23.8 


42.9 


41.5 




38.0 


45.9 


45.5 




Current marital status 


















Currently 


57.7 


59.3 


59.2 


<0.001 


56.2 


63.4 


63.0 


<0.001 


married 


















Marriage 


18.1 


28.7 


27.9 




23.4 


26.4 


26.3 




dissolved 


















Never married 


24.2 


12.0 


12.9 




20.3 


10.2 


10.7 




Age at sex debut, y 


















<15 


21.3 


27.0 


26.6 


<0.001 


25.2 


28.0 


27.8 


0.147 


15+ 


78.7 


73.0 


73.4 




74.8 


72.0 


72.2 




Duration in sex work, y 


















<2 


15.3 


17.3 


17.2 


0.197 


17.4 


16.5 


16.5 


0.011 


2-A 


31.9 


33.3 


33.2 




38.7 


34.7 


34.9 




5-9 


29.1 


26.5 


26.7 




27.7 


26.8 


26.9 




10+ 


23.7 


22.8 


22.9 




16.3 


22.0 


21.7 




Place of solicitation 


















Home 


12.4 


21.5 


20.8 


<0.001 


8.7 


15.4 


15.1 


<0.001 


Brothel/Lodge/ 


43.6 


21.7 


23.3 




29.3 


19.6 


20.1 




Dhaba 


















Public places 


43.2 


53.8 


53.0 




54.6 


58.1 


57.9 




Other 


0.8 


3.0 


2.9 




7.4 


6.9 


6.9 




Experienced violence 


















No 


90.7 


87.8 


88.0 


0.009 


89.6 


92.0 


91.9 


0.039 


Yes 


9.3 


12.2 


12.0 




10.4 


8.0 


8.1 




Has a nonpaying main partner 
















No 


51.0 


36.8 


37.9 


<0.001 


46.0 


37.5 


37.9 


O.001 


Yes 


49.0 


63.2 


62.1 




54.0 


62.5 


62.1 




Consistent condom use with clients 
















No 


21.4 


35.9 


34.8 


<0.001 


15.9 


16.0 


16.0 


0.952 


Yes 


78.6 


64.1 


65.2 




84.1 


84.0 


84.0 




Client volume 


















<10 per week 


53.7 


63.2 


62.8 


<0.001 


52.3 


59.6 


59.2 


O.001 


10+ per week 


46.3 


36.8 


37.2 




47.7 


40.4 


40.8 





*P value for the x test comparing the distribution of characteristics of participants providing and those not providing a biological sample, with 
appropriate degrees of freedom 



The 2 latter program indicators were also highly correlated together 
(Spearman r = 0.647, P = 0.006). The difference in the propor- 
tion of FSWs contacted on a monthly basis between 2006 
and 2008 was not significantly correlated with the changes in 
HIV prevalence (Spearman correlation coefficient = —0.1992, 
P = 0.351) or with the changes in the other program indicators 
(data not shown). 

The left columns of Table 4 show the results of the mul- 
tilevel model without program variables and interaction terms. 
The adjusted odds ratio (AOR) comparing HIV prevalence in 
round 2 with that in round 1 was 0.76, with a narrow 95% con- 
fidence interval (95% CI) of 0.67 to 0.85. A number of demo- 
graphic and behavioral variables were associated with HIV 
(Table 4). Surprisingly, women reporting consistent condom use 



with clients and carrying condoms with them had a significantly 
higher HIV prevalence. The only significant contextual variable 
among those examined was mean age at marriage of girls, with 
an approximately 8% lower HIV prevalence for each year of 
increase in this district-level variable. As expected, baseline 
HIV prevalence was also strongly associated with overall HIV 
prevalence. The random coefficient for the round had a variance 
significantly higher than zero, reflecting the high level of hetero- 
geneity in the changes in HIV prevalence across districts. The 
variance of the random intercept, although much lower than in the 
null model, was also significantly different from zero (Table 4). 

The right columns of Table 4 show that there was a sig- 
nificant negative interaction between the round and changes in 
program coverage between 2006 and 2008. The interaction terms 
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TABLE 2. Comparison of Selected Individual Characteristics of 
Female Sex Workers Between Survey Rounds 1 and 2 



Round 1 Round 2 
Characteristics (n = 10,806) (n = 10,901) P* 

Current age, y <0.001 

15-24 20.8 17.8 

25-34 45.8 44.3 

35+ 33.3 38.0 

Literacy <0.001 

Illiterate 65.2 60.6 

Literate 34.8 39.4 

Source of income other than <0.001 
sex work 

No 57.1 54.1 

Yes 42.9 45.9 

Current marital status <0.001 

Currently married 59.2 63.4 

Marriage dissolved 28.7 26.4 

Never married 12.1 10.2 

Age at sexual debut, y 0.150 

<15 27.1 28.0 

15+ 72.9 72.0 

Duration in sex work, y 0.063 

<2 17.3 16.5 

2-4 33.3 34.7 

5-9 26.5 26.8 

10+ 22.9 22.0 

Place of solicitation <0.001 

Home 21.5 15.5 

BrotheVLodge/Dhaba 21.7 19.6 

Public places 53.7 58.1 

Other 3.0 6.9 

Carrying condoms <0.001 

No 56.0 40.1 

Yes 44.0 59.9 

NG/CT t <0.001 

No 92.8 94.1 

Yes 7.2 5.9 

Experienced violence <0.001 

No 87.8 92.0 

Yes 12.2 8.0 

Has a nonpaying main partner 0.330 

No " 63.1 62.5 

Yes 36.9 37.5 

Syphilis serological test result <0.001 

Negative ' 88.3 92.9 

Positive 11.7 7.1 

Consistent condom use with <0.001 
clients 

No 35.8 16.0 

Yes 64.2 84.0 



*P value for the \ 2 test comparing the distribution of characteristics 
of participants between the study rounds. 

f Infection by either Neisseria gonorrhoeae or Chlamydia trachomatis. 

between the round and changes in other program variables were 
not statistically significant and thus not kept in the model. Both 
baseline HIV prevalence and having a regular nonpaying partner 
also had significant multiplicative interactions with survey round. 
The only variable not involved in the interaction terms that did not 
remain significant compared with the previous model was mean 
age at marriage of girls at the district level. The variance of the ran- 
dom coefficient for the round, though significantly different from 
zero, was lower than that in the previous model, whereas the variance 
of the random intercept was now close to null in this model (Table 4). 

Table 5 shows the decreases in HIV prevalence for the 
various categories of variables shown to have an interaction with 



the study round from the results presented in Table 4. The decrease 
in HIV prevalence was not as strong among women with a main 
nonpaying partner compared with those without, although the 
latter had a higher overall HIV prevalence. The decrease in HIV 
prevalence was not significant in districts with very low baseline 
HIV prevalence and significant for higher values starting at 
about the first quartile of this variable. Finally, there was a 
stronger decline in HIV prevalence in districts with a larger 
increase in program coverage between 2006 and 2008 (Table 5). 
The decrease in HIV prevalence was not significant when 
coverage was higher in 2006 than in 2008. 

DISCUSSION 

This multilevel analysis shows significant declines in HIV 
prevalence among FSWs in southern India between 2005 and 
2009 that are strongly associated with increases in Avahan pro- 
gram coverage, thus suggesting a positive impact of the FSW 
intervention component of this intervention where implemented 
in a way ensuring sufficient increase in program coverage up 
until 2008. Indeed, larger increases in program coverage at the 
district level led to larger declines in HIV prevalence, with no 
significant decrease in districts occurring without an increase 
in coverage, and a reduction of approximately 50% in HIV prev- 
alence in the districts with the highest increases in coverage. 

This analysis controlled for a number of individual and 
contextual factors that could have potentially confounded this as- 
sociation, as well as for baseline HIV prevalence at the district 
level and its interaction with the round, thus taking into account 
the fact that decreases in prevalence were more likely to occur in 
districts with higher baseline values. The use of a random coef- 
ficient for the round took into account the wide variability in the 
changes in HIV prevalence between districts. 

It was difficult to fully address the impact of the program 
indicator on condom requirement met because of the high corre- 
lation between its 2006 to 2008 district-level variation and that 
of program coverage. However, the borderline significant corre- 
lation between increases in condom requirement met and de- 
creases in HIV prevalence at the district-level suggests that this 
indicator is also part of the mechanism leading to impact. Finally, 
it is possible that the actual coverage of the program is more 
important in terms of impact than the frequency of the contacts 
with the target population, as long as some (unknown) threshold 
is reached for the latter. 

The associations that we found between HIV prevalence 
and the individual variables considered in the analysis are largely 
similar to what was found in an analysis of round 1 data. 14 That 
condom use and condom carriage were both associated with a 
higher HIV prevalence may be due to reverse causation, with FSWs 
having acquired HIV before being reached by the interventions 
and starting to use condoms once they learned more about HIV, 
especially when they discovered their serological status. A sim- 
ilar phenomenon was observed in studies of the general popu- 
lation in southern India. 13 ' 16 Among the contextual district-level 
variables considered, only mean age at marriage of girls was 
significantly associated with HIV in the basic model. Interest- 
ingly, this factor was also associated with HIV prevalence among 
pregnant women attending antenatal care clinics in the same 
IBBA districts. 17 However, in the model with program variables, 
mean age at marriage was not a significant predictor of HIV 
Careful examination of the data showed that it was the introduction 
of the interaction between baseline HTV prevalence and the round 
that was mainly confounding this association (data not shown). 

The results of this study constitute a significant addition 
to the overall picture of the Avahan impact evaluation. Whereas 



Sexually Transmitted Diseases • Volume 41, Number 6, June 2014 383 
Copyright © 2014 by the American Sexually Transmitted Diseases Association. Unauthorized reproduction of this article is prohibited. 



Alary et at 



TABLE 3. Comparison of HIV Prevalence (in percent), Program Coverage* and Condom Requirement Met 1 Between 2 Cross-Sectional 
Surveys Carried Out in 2005 to 2006 (Round 1 ) and 2008 to 2009 (Round 2) Among FSWs in 24 Districts of Four Southern States of India 

HIV Prevalence Program Coverage Condom Requirement Met 

Round Round Round 2 Minus Round Round Round 2 Minus Round Round Round 2 Minus 

State/District 1 2 P Round 112 Round 1 1 2 Round 1 



Andhra Pradesh 


16.3 


11.3 


<0.001 


-5.0 


28.1 


81.3 


+53.2 


19.7 


97.4 


+77.7 


Chitoor 


8.7 


12.1 


0.123 


+3.4 


29.8 


96.1 


+66.3 


1.3 


66.5 


+65.2 


East Godavari 


24.4 


21.7 


0.356 


-2.7 


21.5 


52.1 


+30.6 


27.1 


131.6 


+104.5 


Guntur 


20.5 


13.8 


0.012 


-6.7 


31.9 


67.5 


+35.6 


17.8 


132.3 


+114.5 


Hyderabad 


14.8 


13.7 


0.665 


-1.1 


26.2 


90.9 


+64.7 


4.0 


46.6 


+42.6 


Karimnagar 


22.8 


7.0 


<0.001 


-15.8 


30.8 


100.2 


+69.4 


21.4 


71.6 


+50.2 


Prakasham 


10.2 


8.8 


0.519 


-1.4 


26.6 


83.2 


+56.6 


49.2 


133.7 


+84.5 


Visakhapatnam 


14.8 


9.8 


0.027 


-5.0 


26.4 


51.3 


+24.9 


41.0 


136.6 


+95.6 


Warangal 


13.9 


3.7 


<0.001 


-10.2 


20.2 


91.0 


+70.8 


10.2 


41.0 


+30.8 


Maharashtra 


22.5 


22.2 


0.793 


-0.3 


80.8 


80.1 


-0.7 


20.2 


43.8 


+23.6 


Kolhapur 


33.0 


27.4 


0.292 


-5.6 


45.3 


114.1 


+68.8 


28.6 


103.0 


+74.4 


Mumbai 


21.5 


23.0 


0.376 


1.5 


92.9 


87.0 


-5.9 


24.6 


46.9 


+22.3 


Parbhani 


15.3 


14.9 


0.883 


-0.4 


63.2 


84.8 


+21.6 


50.6 


80.0 


+29.4 


Pune 


35.1 


23.2 


O.001 


-11.9 


5.2 


22.3 


+17.1 


3.7 


32.8 


+29.1 


Thane 


12.2 


20.7 


<0.001 


+8.5 


106.0 


97.1 


-8.9 


16.4 


29.6 


+13.2 


Yevatmal 


37.3 


26.8 


0.047 


-10.5 


26.8 


53.2 


+26.4 


40.2 


160.0 


+119.8 


Tamil Nadu 


8.1 


7.2 


0.285 


-0.9 


78.5 


94.2 


+15.7 


88.8 


71.7 


-17.1 


Chennai 


3.2 


2.0 


0.303 


-1.2 


86.5 


106.3 


+19.8 


67.8 


89.9 


+22.1 


Coimbatore 


6.1 


6.8 


0.705 


0.7 


84.8 


93.7 


+8.9 


99.3 


67.2 


-32.1 


Dharmapuri 


12.8 


7.4 


0.011 


-5.4 


106.4 


117.3 


+10.9 


102.3 


82.8 


-19.5 


Madurai 


5.7 


8.3 


0.149 


2.6 


72.4 


88.0 


+15.6 


94.1 


75.2 


-18.9 


Salem 


12.7 


11.3 


0.544 


-1.4 


61.7 


79.6 


+17.9 


86.0 


47.7 


-38.3 


Kamataka 


18.0 


13.0 


O.001 


-5.0 


34.6 


99.7 


+65.1 


30.6 


87.7 


+57.1 


Bangalore 


11.3 


7.7 


0.022 


-3.6 


19.4 


125.6 


+106.2 


35.2 


102.7 


+67.5 


Belgaum 


34.2 


27.5 


0.048 


-6.7 


17.5 


57.5 


+40.0 


13.2 


77.2 


+64.0 


Bellary 


15.6 


14.4 


0.614 


-1.2 


94.2 


111.4 


+17.2 


69.2 


83.8 


+14.6 


Shimoga 


10.5 


9.4 


0.586 


-1.1 


77.5 


120.1 


+42.6 


23.0 


80.7 


+57.7 


Mysore 


24.2 


11.1 


<0.001 


-13.1 


64.7 


97.4 


+32.7 


16.1 


53.8 


+37.7 


Total 


17.0 


14.2 


O.001 


-2.8 


55.1 


87.6 


+32.5 


34.4 


73.6 


+39.2 



*Percentage of the FSW population contacted at least once by the program in the previous year. This was computed from the number of FSWs 
contacted at least once by the program in the previous year divided by the estimation of the size of the FSW population on the same year. 

^Percentage of the condoms required by FSWs for all sex acts with clients in a given year. This was computed based on the number of condoms 
distributed to FSWs by the NGOs in a given year as the numerator, whereas the denominator was the expected number of condoms needed by FSWs 
in the district, based on the estimated number of FSWs in the district and the mean number of client-contacts estimated from the cross-sectional 
survey data. 



it has been shown that the implementation of the program was 
highly successful, 12 that favorable time trends in intermediate 
outcomes (e.g., condom use and prevalence of curable sexually 
transmitted infections) occurred after implementation, 18 " 20 that 
condom use was positively associated 21 and sexually transmit- 
ted infection prevalence negatively associated 22 with exposure 
to the intervention, and that the program seemingly had an impact 
on HIV prevalence at the general population level, 8,23 the evidence 
for an impact on HIV prevalence among FSWs has so far been 
limited to time trends observed 19 - 20 and mathematical modeling 
studies 5 " 7 based on these trends. This is the first study showing 
a direct relationship between program exposure assessed through 
routine program monitoring and decrease in HIV prevalence in 
one of the vulnerable populations targeted by the intervention. 

However, there are a number of limitations to this study. 
First, there was no randomization of the intervention and the 
assessment of the impact relies only on serial observational data. 
The decision not to use a randomized design for this evaluation 
was taken right at the start for reasons that have been summa- 
rized in previous publications. 1 ' 3 The limitations of community- 
randomized trials for "real-world" evaluations have since that 
time been acknowledged more widely. 2 Furthermore, we were 
able to demonstrate an association between the level of inter- 
vention coverage and decreases in HIV prevalence at the district 



level, thus reaching a level of plausibility in the attribution of 
the impact on HIV prevalence to Avahan as per the framework 
proposed by Habicht et al. 24 

Second, we do not know the characteristics of the women 
who declined participation in the surveys and the proportion of 
refusals increased between the 2 rounds. It is impossible to assess 
the impact of refusals on the observed trends because no infor- 
mation is available on the women who did not participate. In ad- 
dition, there were significant differences between women who 
provided biological samples and those who only accepted to be 
interviewed. However, for most variables, these differences were 
of small magnitude, given the very large sample sizes, or went 
in the same direction in both survey rounds. Moreover, the dis- 
tribution of the variables considered in Table 1 among all par- 
ticipants was almost identical to that among those tested for 
HIV, because the proportion of FSWs not providing a sample 
among the participants was low in both rounds. Thus, important 
biases resulting from these differences when comparing HIV 
prevalence in both rounds are unlikely. Similar observations 
are relevant to the comparison of the characteristics of the par- 
ticipants between the rounds (Table 2): apart for the changes 
expected because of the intervention, the absolute differences 
were mostly small, although often statistically significant because 
of the very large sample sizes. 
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TABLE 4. Results of the Multilevel Logistic Regression Models of Individual-, Programmatic-, and District-Level Determinants of HIV 
Prevalence Among FSWs in 24 Districts in South India, 2005 to 2009 



Without Program Variables With Program Variables and 

Characteristics and Interaction Terms Interaction Terms 



Fixed Part of the Model 


AOR* 


95% CI* 


AORt 




95% CI^ 


Intercept 


A 1 O 

U. Iz 


A AT A A/1 

U.Uz— U.04 


A AQ 

U.Uo 


A AAO 

u.uuz 


A AO A A A 

U.Uz— U.4U 


IrSrSA round 












Round 1 


Reference 


Reference 


Reference 


Reference 


Reference 


Round 2 


A 7A 

u. /o 


A A 7 ACS 

u.o /— U.oj 


1 1 o 

i . iy 


A OIA 

U.zjU 


A OA 1 SO 

u.yu— i jo 


Current age, y 












<Zj 


Reference 


Reference 


Reference 


Reference 


Reference 


OS 1/1 
zj — j4 


1 1 0 
1 . 1 o 


1 AS 1 11 

1 .Uj— 1 .jj 


1 1 Q 

i . iy 


A OAS 
U.UUj 


1 AS 1 11 

1 .UJ— 1 .jj 


1 s_i_ 
J J+ 


A Ol 


A Q 1 1 AT 


A Ol 

u.yj 


A 1 1 1 

U.J 1 1 


A Q 1 1 A7 

U.o 1—1 .u / 


Literacy 












Illiterate 


Reference 


Reference 


Reference 


Reference 


Reference 


Literate 


A 07 

U.o / 


A OA A OS 

u.ou— u.yj 


A OQ 

U.oo 


A AA/1 

U.UU4 


A OA A OA 

u.oU— u.yo 


Source of income other than sex work 












INO 


Reference 


Reference 


Reference 


Reference 


Reference 


Van 

Yes 


A OA 

u.yu 


A OO A oo 

u.oz— u.yy 


A OQ 

U.oo 


A A 1 A 

U.U1U 


A OA A 07 

u.oU— u.y / 


Current marital status 












Currently married 


Reference 


Reference 


Reference 


Reference 


Reference 


Marriage dissolved 


1 SA 
1 . JU 


1 11 1 AS 

1. J /—I .to 


1 /I O 

1 .40 


^-A AA1 

<U.UU1 


1 IS 1 A1 
1 . J J— 1 .OJ 


Never married 


l i A 

1 . 1 o 


in? i 11 

1 .uz— 1 Jj 


1 1 s 

1 . 1 J 


n nio 

U.OJ 7 


1 .ul — 1 . J 1 


Age at sex debut, y 












<15 


Reference 


Reference 


Reference 


Reference 


Reference 


1 J+ 


A 01 

U.y 1 


A 01 A no 

u.oj— u.yy 


A 01 

u.y i 


A AOO 

u.uzy 


A 01 A OO 

u.oj— u.yy 


Duration in sex work, y 












<z 


Reference 


Reference 


Reference 


Reference 


Reference 


O A 

z — *t 


1 9£ 
1 .zo 


1 in 1 41 
1 . 1 U — 1 -H-J 


1 97 
1 .z / 




111 1 AA 


s o 

j — y 


1/10 

1 .4o 


1 OO 1 AO 

i .zy— i.oy 


1 A O 

i .4y 


^-A AA1 
<U.UU1 


1 1A 1 7A 
l.jU— 1. /U 


1 0_i_ 


1 £.9 
1 .OZ 


1 A(\ 1 88 
1 .4U— 1 .oo 


i £4 

1 .04 


<-A AA1 
^U.UUl 


] A] 1 OQ 
1 .41 — 1 0 7 


Place of solicitation 












Home 


Reference 


Keterence 


Keterence 


Reference 


Keterence 


Brothel/Lodge/Dhaba 


1 TO 

1. /z 


1 A O 1 OO 

i .4y— iyy 


1 71 

i./j 


^-A AA1 

<U.UU1 


1 SA O AA 
l.jU— Z.UU 


Public places 


1 1A. 

J.jO 


1 90 1 SA 
1 .zu — 1 .Jt 


i i£. 

1 .JO 




1 90 1 S4 
1 .zu — 1 .Jt 


Other 


1 A1 
1 .Uj 


A O 1 111 

U.ol— l.jl 


1 AS 

I.Uj 


A 71 1 
U. / 1 1 


A QO 1 11 

U.oz— 1 .jj 


Carrying condom 












INO 


Reference 


Reference 


Reference 


Reference 


Reference 


Yes 


1 OA 

1 .zU 


1 1 A 111 

1 . 1 u — 1 . j 1 


1 1 A 

l.lo 


A AA1 

U.UU1 


1 A7 1 0 7 

l.U/— l.z / 


JNU/L 1 












Van 

Yes 


1 is 


117 ISA 

1 . 1 /—l .JO 


1 1A 
1 .JO 


<-A AA1 
<~U.UU1 


117 1 S7 
1 . 1 /— 1 .J / 


INO 


Reference 


Reference 


Reference 


Reference 


Reference 


Lxpenenced violence 












INO 


Reference 


Reference 


Reference 


Reference 


Reference 


Yes 


1 1C 

1 Jj 


1 1 O 1 SI 

i . iy— i . j j 


1 IS 

i.jj 


■r-A AA1 

<U.UU1 


1 1 O 1 SI 

i . iy— i . j j 


Has a nonpaying main partner 












INO 


i/ii 

1 .4 I 


1 OO 1 S/l 

i .zy— i .j4 


1 oo 


■r-A AA1 

<U.UU1 


1 1 S 1 A S 

1 . 1 J— 1 .4j 


Yes 


Reference 


Reference 


Reference 


Reference 


Reference 


Syphilis serological test result 












Negative 


Reference 


Reference 


Reference 


Reference 


Reference 


Positive 


Z.J J 


O AO 0 AA 

z.uy— z.ou 


O 10 
Z. JZ 


■r-A AA1 

<U.UU1 


O AO o so 

z.uo— z.oy 


Consistent condom use with clients 












No 


Reference 


Reference 


Reference 


Reference 


Reference 


Yes 


1.19 


1.07-1.32 


1.19 


0.001 


1.07-1.31 


Baseline (round 1) HIV prevalence 8 


1.05 


1.04-1.07 


1.05 


<0.001 


1.04-1.07 


Distal variable 












Mean age at marriage of girls 


0.92 


0.85-0.998 


0.93 


0.076 


0.87-1.01 


Program variables from CMIS^' 












FSW coverage in 2006 






1.00 


0.480 


0.99-1.01 


Difference in coverage between 






1.00 


0.201 


0.99-1.01 


2008 and 2006 












Interaction between round and 












Having a nonpaying partner 






1.24 


0.010 


1.05-1.45 


Difference in coverage between 






0.995 


0.006 


0.991-0.998 


2008 and 2006 












Baseline (round 1 ) HIV prevalence 5 






0.98 


0.001 


0.97-0.99 



(Continued on next page) 
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TABLE 4. (Continued) 



Characteristics 



Random Part of the Model 



Without Program Variables 
and Interaction Terms 



Estimated 
Variance 



95% CI of 
Variance 



With Program Variables and 
Interaction Terms 



Estimated 
Variance 



95% CI of 
Variance 



Estimated district-level variance of 

Round 0.037 
Intercept 0.004 



0.017-0.080 
0.000002-0.035 



0.026 
0.001 



0.011-0.059 
Not significantly 
different from 0 



The random intercept variance decreased from 0.4094 (95% CI, 0.226-0.741) in the null model to 0.004 (95% CI, 0.000002-0.035) in the model 
with a random effect for the round, but no program variables or interaction terms. 

*Multilevel logistic regression model excluding program variables and interaction terms; —2 log likelihood ratio test: \ 2 = 97.54, P < 0.001. 
^Multilevel logistic regression model including program variables and interaction terms; —2 log likelihood ratio test: x 2 = 67.26, P < 0.001. 
^Infection by either Neisseria gonorrhoeae or Chlamydia trachomatis. 
§ HIV prevalence % among FSWs in round 1 of IBBA. 

^ORs for every percentage increase in coverage or one percentage increase in difference in coverage. 
CMIS indicates computerized monitoring information system. 



Third, given the potential of antiretroviral therapy (ART) 
to decrease HIV transmission, 25 the increase in ART coverage be- 
tween the study rounds could have been involved in the decrease 
in HIV prevalence we observed, whereas it could not been directly 
assessed because information on ART use by HIV-infected women 
was not collected in the IBBA surveys. This favorable impact of 
ART on HIV prevalence is, however, unlikely in our study. In- 
deed, out of an estimated 2.6 million HIV-infected people in 
India, only 6845 received ART in 2005, whereas coverage in- 
creased to 140,654 people on treatment (approximately 5% of 
all infected people) in 2008 and to 223,223 (approximately 9% 
of all infected people) in 2009. 26 Antiretroviral therapy coverage 
thus remained relatively low until 2008 to 2009, by which time 
IBBA surveys among FSWs were completed in all districts. In 
addition, as ART improves survival, its expansion could have 



first led to an increase in HIV prevalence, despite a reduction in 
new HIV cases, as was recently observed in a highly affected 
community in South Africa. 27 Consequently, if the increase in 
ART coverage between the study rounds had any impact on HIV 
prevalence trends during this period, it would likely have been to 
artificially mask a downward trend. 

Finally, the first survey was carried out after more than 
1 year of intervention in most districts, and the level of intervention 
coverage had already reached relatively high levels at that time; this 
could have led to an underestimation of the decrease of HrV prev- 
alence between the study rounds and of the association between 
increases in program indicators and declines in HIV prevalence. 

In conclusion, this study strongly suggests that the Avahan 
intervention, where appropriately implemented, had a favorable 
impact on HIV prevalence among FSWs in the 4 southern states 



TABLE 5. Changes Observed in HIV Prevalence Among FSWs Between the Study Rounds (2005-2006 and 2008-2009) at Different 
Levels of the Independent Variables Having a Significant Statistical Interaction With the Study Round*, 24 Districts of South India 

Factors Having a Significant P value of the 

Interaction With the Round AORt 95% CI+ P Interaction Term 



Having a main nonpaying partner 

Yes ' 0.84 

No 0.68 

Baseline HIV prevalence* 

Minimum value (3.2%) 1 .00 

1st quartile (11.3%) 0.86 

2nd quartile (14.8%) 0.81 

3rd quartile (21.5%) 0.72 

Maximum value (37.3%) 0.54 

Difference in FSW coverage between 
2008 and 2006 § 

Minimum value (-8.9%) 0.94 

1st quartile (11.0%) 0.85 

2nd quartile (24.9%) 0.78 

3rd quartile (56.6%) 0.66 

Maximum value (lOe.3% 11 ) 0.5 1 



0.012 



0.72-0.97 


0.02 


0.59-0.79 


<0.001 


0.80-1.25 


0.996 


0.74-0.998 


0.047 


0.71-0.92 


<0.001 


0.63-0.81 


<0.001 


0.42-0.69 


<0.001 


0.77-1.16 


0.564 


0.73-0.98 


0.024 


0.70-0.89 


<0.001 


0.57-0.77 


<0.001 


0.37-0.69 


<0.001 



0.003 



0.008 



The results presented in this table differ slightly from those in Table 4 because here, for the sake of simplicity, we used 3 different models with only 
1 interaction variable in each of them (whereas in Table 4, we used 1 model including all the interaction terms); however, the AOR corresponding to 
each of the main variable and its interaction with round was similar in the 3 separate models and the model used for Table 4. 

1 AOR comparing HIV prevalence in round 2 with that in round 1 at different levels of each independent variable interacting with the round (ad- 
justed for all the independent variables listed in Table 4). 

1 95%CI of the AOR. 

5 At the district level. 

^Coverage can be more than 100% because the denominator of this indicator is based on the estimate of the size of the FSW population at any 
moment in the district, whereas the number of women actually covered in the district could be higher than the denominator in places with high turnover 
of FSWs (mean duration of presence of the FSWs in the district of less than 1 year). 
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of India covered by the intervention. The triangulation of our 
results with those from other approaches used in the Avahan 
evaluation 5 " 8 ' 23 suggests a major overall impact of this program 
on the HIV epidemic in southern India. Interventions targeting 
most at-risk populations, such as FSWs, should be scaled up in 
all countries where specific vulnerable populations contribute 
significantly to the spread of HIV 
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